Reducing entanglement with symmetries: application to persistent currents in 

impurity problems 
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We show how canonical transformations can map problems with impurities coupled to non- 
interacting rings onto a similar problem with open boundary conditions. The consequent reduction 
of entanglement, and the fact the density matrix renormalization group (DMRG) is optimally suited 
for open boundary conditions, increases the efficiency of the method exponentially, making it an 
unprecedented tool to study persistent currents. We demonstrate its application to the case of the 
one-channel and two-channel Kondo problems, finding interesting connections between the two. 



I. INTRODUCTION 

The transport and non-equilibrium behavior of 
strongly correlated quantum many-body systems is one 
of the most challenging problems in condensed matter 
physics. Interactions can give rise some complex and 
intriguing phenomena, such as the Kondo effect, with 
counter-intuitive transport properties in systems with 
impurities^ Thanks to advances in nanofabrication, ex- 
perimentalist can routinely manufacture nanostructures 
that resemble artificial atoms -quantum dots- that can 
be manipulated with an extreme degree of control^, re- 
alizing an ideal playground to test transport theories. 
Understanding the physical phenomena arising in these 
systems has a fundamental technological interest, since 
it could lead to the development of the next generation 
of electronic devices. 

The study of persistent currents in mesoscopic systems 
dates back to the 1980's 6 -. A particular problem that 
has attracted a great deal of interest is the behavior of 
persistent currents through quantum-dot devices^. They 
can shed light on the Fermi liquid properties of these 
systems^, and provide a sorely needed tool to calculate 
conductance^ 8 -, and other transport properties^—. In- 
tcrferomctry devices can also serve as probes for frac- 
tional statistics^, and spin-charge separation^. Unfor- 
tunately, the lack of well controlled analytical methods 
that can deal with Kondo physics, especially in a ring ge- 
ometry, can be a source for disagreement between differ- 
ent theoretical treatments ^ 16 ' 18 . Since quantum Monte 
Carloi^ and Numerical Renormalization Group (NRG)^ 
cannot be applied to the case of a magnetic flux threading 
the ring (appearing as a complex phase in the hopping), 
it would be highly desirable to count with a reliable nu- 
merical technique to treat this problem. 

The DMRG method^ could in principle overcome 
these limitations. However, its efficiency with periodic 
boundary conditions (PBC) is undermined by the struc- 
ture of the quantum entanglement in a ring geometry, 
and its application has been limited to the case of spin- 
less fermionsS^, or to small rings. If we consider a parti- 
tion of our system into two disjoint parts, the quantum 
entanglement between the two subdivisions can be quan- 
tified by the von Neumann entropy S. The von Neu- 
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FIG. 1: Cartoons showing the impurity or quantum-dot (QD) 
coupled to a ring (a) and the equivalent systems with open 
boundary conditions: (b) with the impurity in the center, and 
(c) at the edge. 



mann entanglement entropy determines the number of 
states m necessary to efficiently represent the ground- 
state of the system using a matrix product state (MPS)22, 
m ~ exp (S). Since the entanglement entropy for a sys- 
tem with periodic boundary conditions is twice the one 
for open boundary conditions (OBC)^, a much larger ba- 
sis is needed to simulate it. Besides this fundamental fact, 
the structure of the MPS used by the DMRG algorithm 
as an ansatz to approximate the ground state does not 
properly account for the entanglement introduced by the 
closed ends of the ring2^. This is a problem in transport 
calculations, since the current through the ring scales in- 
versely with the length and the coupling strength in the 
problem, and high precision is required to study large 
systems. Ingenious tricks have been introduced^, but 
their application to general situations is limited. 

Here we revisit a well known canonical transformation 
that was originally introduced in the context of quan- 
tum impurity problems^ 6 -, and referred to as a "folding" 
transformation^, mapping a Hamiltonian with periodic 
boundary conditions onto an equivalent model with open 
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FIG. 2: (a) Error in the energy using a ring with PBC, and 
the equivalent system with OBC, relative to results for OBC 
using m — 1000 states. (b) Entanglement entropy for different 
for different cuts along the chain using both systems. The 
impurity is in the center. All results are for a system with 
L = 20 (L ring = 41), J K = 1, <j> = ty/2, at half-filling. 



boundary conditions. This has two important implica- 
tions: (i) it reduces the entanglement by half, allowing 
for a more efficient representation in terms of MPS, and 
(ii) the OBC enable the DMRG algorithm to simulate 
this problem optimally. 
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FIG. 3: Persistent current as a function of the magnetic 
phase, for (a)(b) one-channel, and (c)(d) two channel Kondo 
rings. All simulations are at half-filling. We show results for 
different values of the Kondo interaction Jk and system sizes. 



II. MODELS AND METHOD 

A. One-channel Kondo problem 

To illustrate the method, we consider, as an illustra- 
tion, the case of a side-coupled Kondo impurity con- 
nected to a one-dimensional non-interacting chain, as pic- 
tured in Fig. 1(a). Let us split the chain into a left and 
right halves, and assume for the moment that the leads 
have both a finite length L. The Hamiltonian is written 
as: 

H -^lcads ~t~ -^boundary 

L-l 

Pleads = ~t ^ S ( C L> C AJ+1, CT + h - C ') (!) 
A=0,1 j = l,<r 

^boundary = ~t ^ ^ ( C \,l.q C 0^ + n - C -) 

A = 0.1 cr 

Hdot = JrS- so — hS z 

where c\j !<T is the electron annihilation operator acting 
on site j of lead A (where the values 0,1 correspond to 
left and right leads), with spin a = {t,l}, and S is the 
spin operator acting on the impurity/dot. The spin sq 
on site "0" connects the chain to the impurity via an an- 
tifcrromagnctic exchange parametrized by Jr. We have 
also included a Zeeman field h acting on the impurity 



spin. We now introduce a symmetric (+) and antisym- 
metric (-) combination of operators acting on the left and 
right leads. This is nothing else but an application of the 
reflection symmetry, yielding new even(+) and odd(-) op- 
erators: 

c ±J > = -^(c 0j > ± c h j !<T ). (2) 

This is a simple change of basis, with the new operators 
obeying fcrmionic anti-commutation rules, that yields a 
curious and convenient identity. The hopping term in 
^boundary becomes: 

boundary = ( C +4,<x C °^ + h ' C ) ' ( 3 ) 

(7 

while the hopping term into the (-) leads has cancelled 
identically. As a consequence, the new Hamiltonian will 
consist of an impurity coupled to a single lead (+), and 
a second decoupled lead (-). This means that we can 
just solve the impurity-lead piece of the system indepen- 
dently, leading to a reduction of the Hilbert space by a 
power of one-half. This transformation has been used in 
NRG for decades^. However, to the best of our knowl- 
edge, its application to systems with PBC has been cu- 
riously ignored, maybe because the NRG method works 
for the thermodynamic limits. 

To illustrate the application of the scheme to a case 
with PBC we studied the behavior of the persistent cur- 
rent at half-filling, as a function of the coupling Jk- Let 
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us introduce a hopping term connecting the last sites of 
the two leads, labelled as "L" . We can easily see that the 
hopping term reduces to: 

Apse = -t J~] (n+,L,<r - n-,L,o-) , (4) 

cr 

where n± = c\_c± is the density operator. This term 
is just a boundary chemical potential with opposite signs 
for the (+) and (-) leads. But the leads remain decoupled, 
meaning that we can still solve the problem of a single 
lead with OBC! 

Let us consider now the case of a magnetic flux thread- 
ing the ring. This problem was extensively studied in 
RefQl 

using DMRG in systems with up to L r ing — 35 
sites in an enormous computational effort using PBC. 
The flux (j) is introduced by adding a complex phase in 
the hopping matrix element t — > iexp (icj) / L r i ng ) , where 
the total length of the ring is L r i n g = 2L + 1 in our 
notation. By performing a gauge transformation on the 
fermionic operators, we have the freedom to move the 
phase to any link along the ring. In particular, we are 
going to move it to the one connecting the site "0" to the 
leads in -//boundary in (1)- The term containing the phase 
will now read: 

^boundary = ( e ^4,l,c7 C 0.<r + C [l,/0,» + n - C -) • ( 5 ) 

cr 

The flux in the hopping term now introduces a complica- 
tion: the (-) lead will no longer decouple and, as a conse- 
quence, the new, transformed Hamiltonian will have an 
additional term: 

^boundary = ~ V^* £ (e^ COS (^/2)c^ 1|ff Co )<r 

(7 

+ ie^sinO^cl ^o^ + hx.) . (6) 

This hopping term now couples to both (+) and (-) chan- 
nels, because we have broken reflection symmetry, but 
the system still has OBC, as depicted in Fig. 1(b). 

Fig Ell a) shows the error in the energy a function of 
the number of states for a ring with J% = 1,0 = tt/2, 
L = 20 (L r i n g = 41), using PBC in real space, and the 
equivalent system with OBC in the transformed basis. As 
a reference we have used the energy of the system with 
OBC and keeping m = 1000 DMRG states. We clearly 
see that we can achieve better accuracy with OBC, using 
a fraction of the number of states. This can be explained 
by looking at the behavior of the entanglement entropy, 
shown in Figl^b), where the entanglement entropy of the 
transformed system with OBC is smaller by a factor ~ 2. 
In this plot, we show values of S for different cuts along 
the system, with site "0" and the impurity situated both 
at the center. 

Making use of the gauge transformation, we can choose 
to put the phase factor at the link connecting both leads. 
In that case, the (-) channel would decouple from the 
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FIG. 4: Cartoons for the 2-channel Kondo problem showing 
(a) the impurity coupled to two rings, and (b) the equivalent 
system with open boundary conditions. 

site "0", recovering the expression but the boundary 
term would have to be corrected: 

i/pBC = ~t (cos ((j))n +tL ^ - cos (cf))n-. L . a + 

cr 

+ ism(<f>)c!_ tLt(r c +! L,a + h.c.) . (7) 

This scenario is depicted in Fig. 1(c). 

The current can be obtained by differentiating the en- 
ergy as a function of the flux J = —dE/dcfi. For nu- 
merical convenience, we choose to calculate the expres- 
sion J = — it(cQ 1 a co ~ CqCq,i,u)-i since it can be di- 
rectly obtained from the ground-state. Depending on 
the gauge choice, we measure it either at the link con- 
nect site "0" , or at the one connecting the two leads. In 
FigEJa) we show results for a one-channel Kondo ring 
with L r i n g =81. These results were obtained effortlessly 
using m — 600 states, although the error is small enough 
with m = 200 for intermediate values of Jk- The profile 
of the persistent current can be deduced from symmetry 
consideration a 16 ' 27 . Using reflection symmetry around 
the dot (coj,er — > ci JiCr ), we find that J(4>) = J(—<f>). 
This indicates that it is sufficient to estimate j in the 
interval < <f> < n. Furthermore, for the present case 
with TV = Lring odd, an electron-hole transformation on 
the fermions (c\ t j t(T c \ j &) an< 4 a lso on the impurity, 
maps H(—(f>) onto H(<j> + it) and J(—<fi) = —J{4> + tt). 
Combining both results we obtain that </(</>) = J(0 + 7r), 
and J(n/2) = 0. This means that the persistent current 
will have a periodicity in n, instead of 2-7T. 

With little additional complication, the transformation 
can be applied to a ring with an even number of sites. For 
comparison, we show results for a system with L r i ng = 80 
in FigGHb). Using similar symmetry arguments, one can 
show that the periodicity in this case should be 2ir. In 
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FIG. 5: Current per spin for (a) the one-channel Kondo 
problem with a Zeeman field, and (b) the two-channel Kondo 
problem. 



all cases the current is rapidly suppressed with increasing 
Jk- For large Jk, the size of the "Kondo cloud" will 
be confined to the site "0" in direct contact with the 
impurity. This local tightly bound singlet will suppress 
conduction through the ring. 



and L r i n g = 40 respectively, meaning that the total 
lengths of the systems are L to t = 83 and L to t = 81 sites. 
Interestingly, the current has the same behavior as in the 
one-channel problem with L r i ng odd, even though in this 
case, since we have two rings, the total number conduc- 
tion electrons is always even. For L r i ng = 41, one would 
expect a periodicity in 2tt instead of it. By looking at the 
currents for each spin sector in FigJSJb), we observe that 
they both have a periodicity of 2ir, but they are shifted 
by a phase 7r. Thus, the total current will have a period 
7r. This resembles the situation encountered in a single 
side-coupled quantum dot with an applied magnetic field, 
that was proposed as an efficient spin filter in Refs l30ll3l1 . 
A simple intuitive explanation would be to assume that 
the screening effect of the second ring on the impurity is 
seen by the first ring as an effective magnetic field act- 
ing on the impurity. To prove this picture, we applied a 
small Zeeman field to the one-channel problem, and we 
show the results in FigJSJa). We clearly observe the same 
behavior as in the two-channel case, with curves that are 
qualitatively indistinguishable. 



III. SUMMARY AND CONCLUSIONS 



B. Two-channel Kondo problem 

We now generalize the above considerations to the case 
of a two-channel Kondo model^. In this case, the impu- 
rity is also coupled to a second channel, represented by a 
second ring in FigUJa). We can perform a similar canon- 
ical transformation in both rings independently. We as- 
sume that a magnetic flux </> is threading both rings, and 
we put the phase terms on the connecting link at the end 
of each lead as in Eq.©. This open both rings symmet- 
rically about the impurity site (see Fig0£b)), yielding 
an equivalent one-dimensional system with OBC, that 
can easily and efficiently be simulated with the DMRG 
method. The remarkable aspect in this case is that we 
have eliminated two closed rings from the problem. 

In this case, we study the problem in the sector with 
total S z = 1/2, and even total number of conduction 
electrons. Unlike the one-channel case, in the present 
case the impurity will be overscreened, with each channel 
trying to form a singlet with the impurity. Thus, the 
impurity is expected to be less efficient at suppressing 
the current through the leads. 

In Figs[3{c) and (d) we show results for the persistent 
current through an impurity coupled to two channels. 
The length of each lead, for each channel is L r i ng = 41 



To summarize, we have applied a "folding" canoni- 
cal transformation to lattice models of impurities cou- 
pled to rings, mapping the problems onto equivalent sys- 
tems with open boundary conditions. As a remarkable 
counter- intuitive result, we find that a system with open 
ends can realize persistent current when transformed to 
the original real-space basis. As a consequence, entan- 
glement is reduced by a factor 2, and makes the prob- 
lem suitable for efficient DMRG simulations of large sys- 
tems with little or moderate effort, as shown here for the 
cases of the one-channel and two-channel Kondo prob- 
lems. The folding transformation cannot be applied to 
disordered rings , or to problems with bulk interactions, 
since many-body terms would introduce long-range in- 
teractions between (+) and (-) leads. However, the sig- 
nificant entanglement reduction opens the doors to finite 
temperature calculations^, and the study of complex im- 
purity problems. 
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